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ABSTRACT 


Current,  temperature  and  pressure  measurements  from  instruments 
moored  in  the  MODE  area  during  the  years  1971-1975  are  presented.  Record 
lengths  vary  from  2 weeks  to  7 months.  Data  from  individual  instruments 
are  displayed  in  the  form  of  statistical  tables,  spectral  plots  and 
variables  vs.  time  plots.  Current  da  a are  additionally  displayed  as 
progressive  vector  plots.  Composite  plots  of  temperature,  pressure  or 
velocity  are  included  for  groups  of  records  which  could  be  combined  to 
show  a meaningful  spatial  or  temporal  array.  Selected  CTD  data  are 
included  and  displayed  as  potential  temperature  and  salinity  versus 
pressure. 
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PREFACE 


' 


This  volume  is  the  sixteenth  of  a series  of  data  reports  presenting 
moored  current  meter  and  associated  data  collected  by  the  W.H.O.I.  Buoy 
Group. 

Volumes  I through  XV  present  data  from  the  years  1963-1971,  and 
from  several  special  experiments:  the  1970  Pollard  array,  the  1973  IWEX 
array,  the  1973  MODE  array,  and  the  Saint  Croix  mooring  measurements. 

Volume  XVI  presents  MODE  region  data  from  1971  through  1975. 


i 

i 


Volume  * W.H.O.I.*  * Notes 

Number  * Ref.  # * *year  experiment 

************************************************************************** 


I 
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Webster, 

F. , and 
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II 
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Webster, 

F . , and 

N.  P.  Fofonoff 

III 

67-66 

Webster , 

F. , and 

N.  P.  Fofonoff 
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70-40 

Pollard, 
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Tarbell, 

S . , and 

F.  Webster 

1966 

measurements 

VI 
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Tarbell , 
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S.  Tarbell 

1968 
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VIII 

75-7 
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and  S.  Tarbell 

1970 
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IX 

75-68 

Tarbell , 
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3.  Briscoe, 
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1973 
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76-40 

Tarbell , 
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1969a  measurements 

XI 
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Tarbell, 
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D. , and 

S.  Tarbell 

1973 

MODE  Array 

XIII 

77-18 

Tarbell, 

S . , and 

A.  W.  Whitlatch 

1970 

measurements 

XIV 

77-41 

Tarbell , 

S.,  R.  ] 

Payne,  and 

R.  Walden 

1976 

Mooring  592, 

Saint  Croix 

XV 

/77-56 

Tarbell , 

S . , and 
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1971 

measurements 
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Presentation 


This  report  is  shown  in  its  entirety  on  microfiche  with  a summary 
presented  on  printed  pages  for  easier  access.  The  summary  gives  a 
general  description  of  the  MODE  area  moorings,  some  details  of  the 
instrumentation  used  in  acquiring  the  data  and  a description  of  the 
computer  display  programs.  It  also  contains  composite  plots  of  filtered 
current,  temperature  and  pressure  data. 

The  first  page  of  microfiche  contains  the  summary  and  selected 
CTD  plots.  The  CTD  plots  are  shown  only  on  microfiche.  The  seven  rows 
of  the  first  microfiche  page  are  organized  as  follows: 

Row  A Text 
Row  B Text 

Row  C Composite  plots  by  mooring  location 
Row  D Composite  horizontal  array  plots 
Row  E Composite  temperatures  from  T/Ps 
Row  F Composite  pressures  from  T/Ps 

Row  G Plots  of  temperature  and  salinity  from  CTD  lowerings 
The  basic  current  meter  displays  are  organized  on  subsequent 
microfiche  sheets  so  that  the  displays  for  an  individual  instrument  are 
arranged  down  a microfiche  column,  and  similar  plots  for  different 
instruments  are  arranged  across  the  rows. 

The  seven  rows  of  data  displays  on  microfiche  are  as  follows: 

Row  A Mooring  diagrams 
Row  B Data  identifier 

Row  C Instrument  information  and  computed  statistical  quantities 

Row  D Spectral  diagrams 

Row  E Progressive  vector  plots 

Fow  F u and  v plotted  against  time 

Row  G Speed,  direction  (and  temperature  if  available)  plotted  against  time 
On  microfiche  #6,  the  last  3 columns  have  a different  format  due 
to  the  length  of  the  time  series. 

No  data  are  presented  for  the  various  tensiometers  and  inclinometers 
set  on  some  moorings,  nor  for  experimental  depth  or  temperature  recorders. 
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Introduction 


The  data  presented  in  this  rejwrt  are  from  the  MODE  region  here  defined 
to  be  a circular  area  of  about  1000  square  kilometers  centered  on  28°N, 

69°  40'W  (see  Fig.  1).  The  records  span  3^  years  from  1971  to  1975. 

The  first  array  (MODE-0,  Array  1)  consisted  of  two  surface  and  five 
intermediate  moorings  located  around  and  to  the  west  of  MODE-Center,  (28°N, 

69°  40' W)  where  the  bottom  topography  was  smooth.  The  second  array  con- 
sisted of  four  intermediate,  1 surface  and  a deep  mooring,  located  east  of 
MODE-Center  where  the  bottom  topography  was  rough.  The  third  array  was 
more  widely  spaced,  and  consisted  of  four  surface  moorings,  one  with  a wind 
recorder,  and  one  intermediate  mooring.  The  fourth  array  consisted  of  two 
intermediate  moorings  at  MODE-Center  and  MODE-East  (28°N,  68°  40' W) . The 
fifth  array  (Chausse,  Tarbell,  1976)  was  set  as  part  of  the  MODE-I  experi- 
ment. 

The  arrays  set  after  MODE-I  are  referred  to  as  site  moorings.  Their 
locations  were  at  MODE-Center  (28°N,  69°  40'W)  and  MODE-East  (28°N,  68°  40'W) 
and  intermediate  moorings  were  used  throughout. 

Details  of  mooring  locations  and  durations  at  sea  are  presented  in 
Figures  1,  2 and  Table  1. 

The  first  three  MODE-O  arrays  were  deployed  to  find  an  optimum  design 
for  the  MODE-I  experiment  (Heinmiller,  1975) . They  were  in  different  locations, 
of  different  spatial  density  and  of  different  mooring  types.  The  fourth 
array  and  the  post-MODE  site  moorings  were  deployed  to  obtain  long-term 
records  for  study  of  low  frequency  phenomenon  (Schmitz,  1977) . Mooring  542 
was  set  as  part  of  a subsequent  experiment  and  was  also  a continuation  of 
the  moorings  set  at  MODE-Center.  It  is  included  in  this  report  to  show 
the  full  extent  of  the  data  gathered  at  that  location. 
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MODE  0 ARRAY  2 
FEB.  72-  MAY.  72 


MODE  0 ARRAY  1 
OCT.  71  -FEB.  72 


M • # TO 


M*W  T 0#  69*  68  •* 


SITE  MOORINGS 
OCT.  72  —MARCH  ^3 


MODE  0 ARRAY  3 
MAY.  72-OCT.  72 


MODE  1 ARRAY 
MARCH  73- JUNE  73 


SITE  MOORINGS 
JUNE  73-  APRIL  75 


FIGURE  1 MOORING  POSITIONS 


FIGURE  2 
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Table  1 

Summary  of  MODS 

0 and  MODF 

Site  Moorings 

Moor ing 

Type  * 

Location 

Bottom 

Set, 

Retrieval 

Dates 

Number 

S I IW  B 

( lat . N . , 

long.  W) 

Depth 

(m) 

Array  1 

406 

S 

27°59. 8 ' 

70°  0.3' 

5460 

Oct. 

29 ' 71-Feb. 

7 ' 72 

407 

S 

28°  0.4' 

70°20 . 6 ' 

5465 

Oct. 

30'71-Feb. 

9 '72 

408 

I 

27  °49 . 0 ' 

70°  8.8' 

5470 

Oct . 

30 ' 71-Feb. 

9*72 

409 

I 

28°  1.5' 

70°  6.8' 

5465 

Oct. 

30'71-Feb. 

9 '72 

410 

I 

28  °21 . 5 ' 

69°41 . 5 ' 

5460 

Oct. 

31' 71-Feb. 

9 '72 

411 

I 

28°  0.7' 

69°31. 3 ' 

5427 

Oct . 

31 '71-Feb. 

7 ' 72 

412 

I 

28°  0.2' 

69°41 . 5 ' 

5455 

Oct. 

31' 71-Feb. 

7 ' 72 

Array  2 

424 

S 

28°  9.1 

68° 36. 8' 

5254 

Feb. 

6 '72-May 

25'  72 

425 

S 

28°  0.8' 

69°39.8' 

5462 

Feb. 

8 ' 72-Lost 

430 

B 

28°  9.8' 

68°34 . 5 ' 

5221 

Mar . 

18'72-May 

25 ' 72 

431 

I 

28°20. 3 ' 

68°24 . 8 ' 

5370 

Mar. 

19'72-May 

26 '72 

432 

I 

28°10 . 5 ' 

68°23 . 8 ' 

5380 

Mar. 

20  * 72-Aborted 

433 

I 

28° 10 . 5 ' 

68°23.8' 

5380 

Mar . 

20'72-May 

27 ' 72 

434 

I 

28°10 . 2 ' 

68°11 .8 ' 

5275 

Mar. 

20'72-May 

26 ' 72 

435 

I 

27°58. 7 ' 

68°24 . 7 ' 

5280 

Mar. 

21'72-May 

26 ' 72 

Array  3 

451 

S 

28°54 . 4 ' 

69°41 . 0 ' 

5437 

May 

23 ' 72-Nov. 

1 ' 72 

452 

_ 

27°59. 8' 

70°38 . 7 ' 

5452 

May 

23 ' 72-Oct . 

31 '72 

453 

S 

28°10. 1' 

68°38 . 2 ’ 

5261 

May 

25 ' 72-Oct . 

30 ' 72 

4 54 

S 

27  °33 . 7 ' 

69°41 . 8 ' 

5462 

May 

28 '72-Oct. 

31 ' 72 

455 

s 

28°  0.6' 

69°37 . 6' 

5462 

May 

28 '72-Nov. 

4 ' 72 

Site  Mooring 

473 

I 

28°10. 7 ' 

68°36 . 1 ' 

5261 

Oct. 

30'72-Mar. 

11 ' 73 

474 

I 

28°  1.4' 

69°  39 . 4 ' 

5462 

Nov. 

4'72-Mar. 

10'73 

Site  Mooring 

502 

I 

28°  8.9’ 

68°41 . 4 ' 

5255 

June 

26'73-Dec 

13 ' 73 

503 

I 

28°  0.2' 

69°  44 . 4 ' 

5461 

June 

27'73-Dec 

14 ' 73 

521 

I 

28°  9.5' 

68  °41 . 5 ' 

5265 

Dec . 

13 ' 73-Apr 

20 ' 74 

522 

I 

28°  0.5' 

69°44 . 8 ' 

5462 

Dec . 

14 ' 73-Apr 

21 '74 

538 

IW 

28°  2.6' 

69°44 . 8 ' 

5457 

Apr . 

18 ' 74-Jul 

27 ' 74 

540 

I 

28°  8.9' 

68  °40 . 4 ’ 

5265 

Apr . 

20'74-Jul 

29 ' 74 

542 

IW 

28°  1.4' 

69°  38 . 9 ' 

5462 

July 

27 ' 74-Apr 

;6'75 

* S = Surface 

compound  mooring 

I = Intermediate  compound 

mooring 

IW  = Intermediate  all  wire 

mooring 

B = Bottom  : 

mooring 

1-A-ll 

J 

Moorings 


Both  surface  and  intermediate  moorings  were  set  in  the  early  1970s. 
(See  Table  1.)  Growing  evidence  that  current  meter  data  were  contaminated 
by  noise  from  surface  waves  propagating  down  the  mooring  line  led  to 
the  increasing  use  of  intermediate  moorings.  The  compound  intermediate 
mooring  (Heinmiller,  1976)  proved  extremely  reliable  with  only  3 failures 
in  103  moorings.  In  the  MODE-1  experiment  there  were  large  discrepancies 
between  the  intended  and  actual  instrument  depths  (Chausse  and  Tarbell, 

1976) , which  were  due  in  part  to  the  unpredictable  stretch  factor  of  the 
Dacron  line.  In  1974,  the  first  all-wire  intermediate  moorings  (538, 

542)  were  set  in  an  effort  to  maintain  the  instrumentation  at  its  intended 
depth. 

Numerous  minor  improvements  were  made  over  the  years  to  increase 
the  reliability  of  various  mooring  components.  For  instance,  glass  spheres 
were  placed  in  plastic  hard  hats  instead  of  nets  to  reduce  the  risk  of  an 
implosion  severing  the  mooring  line.  Such  changes  enabled  moorings  to  be 
set  with  confidence  for  9 months  or  longer  as  scientific  necessity  dictated. 

For  greater  detail  about  MODE  area  moorings,  see  Heinmiller  (1975)  , 
Heinmiller  (1976)  and  Berteaux  and  Heinmiller  (1973). 
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Current  Meter  Types 

The  current  meters  described  in  this  report  were  either  Vector 
Averaging  Current  Meters  (VACMs) , built  by  American  Machine  and  Foundry 
(AMF)  or  Model  850  current  meters  built  by  Geodyne,  now  a part  of 
Egerton,  Germeshausen  and  Grier  (EG&G) . 

The  Model  850  current  meter  stores  burst  sampled  data  on  magnetic 
tape  cartridges.  Each  Model  850  collected  and  stored  23  (or  24)  data 
cycles  which  were  sampled  every  5.27  seconds.  The  instrument  then  turned 
off  for  the  remainder  of  the  recording  interval,  usually  either  15  or 
30  minutes.  Those  Model  850s  which  were  modified  to  include  temperature 
measurements  accumulated  the  count  from  the  temperature  circuit  from 
one  5.19  second  period  and  stored  it  at  the  beginning  of  each  data  burst. 

The  VACM  stores  variables  on  1/8"  magnetic  tape  cassettes  at  the 

end  of  each  recording  interval,  typically  every  30  minutes.  Values  in  the 
East  and  North  register  are  computed  from  compass,  van'-  and  rotor  values 
each  time  a rotor  magnet  passes  the  sensing  diode.  The  velocity  components 
are  accumulated  over  the  recording  interval  resulting  in  vector  averaged 

velocities.  The  temperature  count  is  summed  from  the  pulses  of  a voltage 

to  frequency  converter,  whose  output  frequency  is  relayed  to  the  thermistor 
resistance  at  its  input  (Payne  et  al.,  1976). 


Temperature/Pressure  Recorder 

An  instrument  to  measure  and  store  temperature,  pressure  and  time 
data  (T/P)  was  developed  in  the  Draper  Laboratory  at  M.I.T.  for  the  MODE-1 
Data  Array,  and  used  extensively  on  the  post-MODE  moorings.  The  instrument 
stores  a data  sample  every  15  seconds  and  records  the  sum  of  64  successive 
data  samples  on  a magnetic  tape  cassette  every  16  minutes  (64  x 15  = 960  seconds 
= 16  minutes) . 

The  time  base  generator,  a crystal  oscillator,  has  an  accuracy  of  ±1 
second  per  day. 

The  temperature  has  an  accuracy  of  about  .01°C  and  a resolution  of 
•001°C  (Wunsch  and  Dahlen,  1974). 

The  pressure  sensor  is  a strain  gauge  with  a manufacturer  specified 
accuracy  of  .03%  of  the  scale  range  used  (Wunsch  and  Dahlen,  1974)  . 

CTD 

The  CTDs  were  constructed  at  W.H.O.I.  by  Neil  Brown  (Brown,  1974; 

Fofonof f , Hayes  and  Millard,  1975;  Brown  and  Morrison,  1978).  The  instru- 
ments measure  conductivity,  temperature  and  pressure.  Plots  of  potential 
temperature  and  salinity  versus  pressure  are  included  on  the  first  page 
of  microfiche. 
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A time  base  for  the  T/Ps,  850s  and  VACMs  is  set  by  a quartz 
crystal  oscillator  with  a manufacturer's  specified  accuracy  of  ±1 
second  per  day. 

To  provide  secondary  time  information  for  current  meter  records, 
artificial  events  are  placed  in  the  rotor  field  both  before  and  after 
the  sea  data.  The  procedure  for  this  is  as  follows: 

For  First  Events: 


Step 

1. 

Turn  the  recording  circuit  on  (not  to  be  turned 

off  until  recovery) . 

Step 

2. 

a)  Block  the  rotor  to  prevent  it  from  turning. 

b)  Allow  the  instrument  to  record  several  data 

cycles  with  the  rotor  immobilized. 

Step 

3.  (t) 

Ten  seconds  after  a recording  interval  starts,  spin 

the  rotor  as  fast  as  possible. 

Step 

4.  (t) 

In  a successive  record  spin  the  rotor  again. 

Step 

5. 

Block  rotor  again  until  just  before  launch. 

Step 

6.  (t) 

Free  the  rotor  at  launch. 

Last  Events: 

Step 

1.  (t) 

After  recovery,  block  rotor  for  several  recording 

intervals. 

Step 

2.  (t) 

Ten  seconds  after  a recording  interval  starts,  spin 

the  rotor. 

Step 

3.  (t) 

Spin  the  rotor  again  in  a subsequent  record. 

Step 

4. 

Block  the  rotor  for  several  records  before  shutting 

off  the  record  circuit. 

The  steps  that  are  marked  with  a (t)  are  the  most  easily  seen  in  the 
data  and  times  should  be  noted  using  a radio  time  signal  e.g.,  WWV  or  CHU. 

Under  normal  operating  conditions  the  elapsed  time  measured  by  the 
clock  word  matches  the  real  elapsed  time. 

In  this  repcrt  time  is  read  as  year-month-day  hour .minute. second. 
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Current  Meter  Data  Processing 

Bit  strings  from  magnetic  tapes  for  Model  850s  (1/4"  2-track 
cartridges)  and  VACMs  (1/8"  4-track  cassettes)  were  transcribed  onto 
9-track  computer  compatible  tape  at  W.H.O.I.  The  data  were  then  con- 
verted to  scientific  units  (decoded)  and  stored  on  magnetic  tape  in 
Maltais  format  (Maltais,  1969) . 

Editing  the  data  included  selecting  times  for  the  beginning 
and  ending  of  each  record,  adjusting  the  nominal  depth  of  some  records 
to  agree  with  information  supplied  by  the  T/Ps,  applying  corrections 
to  temperature  indicated  by  post  cruise  thermistor  calibrations,  com- 
puting vector  averaged  components  for  the  data  from  burst  sampled 
Model  850  instruments,  and  interpolating  through  gaps  in  the  data 
caused  by  the  removal  of  erroneous  records. 

The  result  is  an  evenly  spaced  time  series  which  is  used  as  a 
basic  input  for  all  further  processing. 

The  current  data  were  classified  into  three  categories  (good, 
questionable  and  bad) . "Good"  is  defined  to  mean  that  all  sensors 
worked  properly.  "Questionable"  means  that  the  data  has  a known  problem 
which  might  not  affect  the  low  frequency  data,  e.g.,  5032.  "Bad" 
means  that  at  least  one  of  the  sensors  did  not  work  at  all.  The  in- 
dividual data  series  presented  in  this  report  contain  mostly  "Good" 
and  some  "Questionable"  data.  Questionable  data  is  indicated  on  the 
Instrument  page  (row  C on  the  microfiche  pages  presenting  data) . 

Table  2 indicates  which  variables  are  presented  with  reference 
to  mooring  numbers  and  approximate  depths. 


Data  Identifiers 


To  insure  that  each  data  series  has  a unique  identifier,  the 
following  guide  lines  are  used;  538,10B900 

538  - The  first  three  digits  are  the  mooring  number 

,10  - The  10th  instrument  down  the  mooring  line.  A one  digit 

mooring  position  is  not  preceded  by  a comma  (5385D900) . 

B - The  position  of  the  letter  in  the  alphabet  indicates  the 

amount  of  editing  that  has  been  done.  The  symbol  $ 
means  no  editing  has  been  done. 

900  - Vector  averaging  interval  (in  seconds) . 1H  would  mean 

1 hour  averages,  ID  for  1 day  averages. 

1DGAU24  - A 1 day  averaged  series  from  Gaussian  filtered  data, 
the  filter  having  a half  width  of  24  hours. 

Table  3,  which  has  the  same  format  as  Table  2 indicates  the  position 
number,  depth  and  type  of  each  instrument  which  has  data  in  this  report. 


Pre  MODE  Data  Series  - MODE  0 

0 

1 

406 

407 

408 

Array  1 
409  410 

Array  2 

411  412  424  430  431  433  434 
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l 

451 

452 

Array 

453 

3 

454 

V 

4 55 

500 

V 

V 

VT 

VT 

VT 

VT 

800 

V 

1500 

VT 

VT 

V 

VT  V 

V VT  VT 

VT 

VT 

1600 

V 
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V 

V 

V V 

VT  VT  V V V 

T 

PT 

T 

4200 

V 

V 

V 

V 

V 

5000+  VT  V 

Post  MODE  Data  Series  - Site  Moorings 


473 

474 

502 

503 

521 

522 

538 

540 

542 

500 

V 

V 

V 

V 

V 

V 

VT 

V 

VT 

SOO 

TP 

TP 

TP 

TP 

800 

TP 

TP 

TP 

TP 

1000 

TP 

TP 

TP 

TP 

1500 

V 

V 

V 

V 

V 

V 

VT 

V 

VT 

2000 

TP 

TP 

TP 

TP 

2500 

TP 

TP 

TP 

TP 

3500 

TP 

TP 

TP 

TP 

4000 

V 

V 

V 

V 

V 

VT 

V 

V 

TP 

TP 

TP 

TP 

4500 

TP 

TP 

TP 

TP 

5000+ 

TP 

TP 

TP 

TABLE  2 VARIABLES  PRESENTED  V=VEL0CITIES 

T-TEMPERATURES 

P=PRESSURE 
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Meters 

0 

500 

800 

1500 

1600 

4000 

4200 

5000+ 

Water 

Depth 

(m) 

500 

600 

800 

1000 

1500 

2000 

2500 

3500 

4000 

4500 

5000+ 
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Instrument  Depth  and  Data  Identifier  Summary 


Pre 

MODE 

- MODE  0 

Array  1 

l 

Array  2 

1 

Array  3 

406 

407 

408  409  410  411 

412  424 

430 

431  433  434 

435  451 

452 

453  454 

0 

(1) 

514 

514 

515 

561 

514 

(3) 

(1) 

(2) 

(2) 

(2) 

816 

(4) 


1518 

(5) 

1516 

(2) 

1503 

(1) 

1522 

(1) 

1504 

(1) 

1476 

(I) 

1502 

(1) 

1519 

(2) 

1516 

(3) 

1620 

(6) 

4003 

(7) 

4001 

(3) 

4028 

(2) 

4008 

(2) 

4074 

(3) 

3975 

(1) 

3955 

(2) 

3990 

(2) 

3970 

(2) 

3968 

(4) 

3970 

-5- 

4075 

(6) 

455 


514 

(2) 


1516 

(3) 


4202 

(8) 


4191  4181  4207  4208 

(3)  (6)  (4)  (4) 

5131  5119 
(4)  (2) 


5460  5465  5470  5465  5460  5427  5455  5254  5221  5370  5380  5275  5280  5437  5452  5261  5462  5462 


Post  MODE  - Site  Moorings 


473 

474 

502 

503 

521 

522 

538 

540 

542 

370 

583 

516 

494 

508 

491 

511 

509 

495 

(2) 

(2) 

(1) 

(1) 

(1) 

(1) 

(1) 

(1) 

(1) 

605 

592 

612 

610 

-2- 

-2- 

-2- 

-2- 

802 

809 

807 

796 

-3- 

-3- 

-3- 

-2- 

999 

1006 

1004 

996 

-4- 

-4- 

-4- 

-3- 

1385 

1595 

1524 

1501 

1506 

1497 

1500 

1511 

1499 

(4) 

(3) 

(2) 

(2) 

(5) 

(5) 

(5) 

(5) 

(4) 

2007 

1998 

2001 

2012 

-6- 

-6- 

-6- 

-6- 

2509 

2495 

2508 

2512 

-7- 

-7- 

-7- 

-7- 

3510 

3497 

3509 

3499 

-8- 

-8- 

-8- 

-6- 

3900 

4105 

4015 

4011 

3998 

3998 

4008 

4000 

(5) 

(4) 

(3) 

(9) 

(9) 

(9) 

<9/ 

(7) 

4014 

4001 

4001 

4003 

-10- 

-10- 

-10- 

-8- 

4414 

4402 

4394 

4409 

-11- 

-11- 

-11- 

-11- 

5166 

5361 

5355 

-12 

-12- 

-12- 

5261 

5462 

5255 

5461 

5265 

5462 

5457 

5265 

5462 

TABLE  3 DEPTHS  OP  INDIVIDUAL  INSTRUMENTS  IN 
METERS 

INSTRUMENT  POSITION  NUMBER  AND 
INSTRUMENT  TYPE;  (N)  for  CM,  -N-  for 
T/P 
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Computer  Programs 

The  following  programs  have  been  used  to  generate  the  data  displays 
in  this  report: 

Statistics  (STATS) 

Standard  statistical  parameters  are  calculated  for  data  in  the  time 
range  given  at  the  bottom  of  the  table.  Given  n speed  and  direction 
or  temperature  values  in  a sample,  we  define  E^  = sinCh,  cosGh  , 

then  for  A = E,  N,  S,  and  T, 


l,  A = - ) A. 

n Z_,  1 


variance 


' A2  - A2 

A n L-.  i 


standard  error  of  the  mean  = — 


standard  deviation  = 0_ 


skewness 


n 

= ^ A? 

n • ! 1 
A L 1=1 


— V A2  + 2A 
n i 

i=l 


1 | 1 \ 4 4A  \ j §A2V  a2 

kurtosis  = — ~ L Ai  “ T L Ai  " i 


Note  that  the  kurtosis  becomes  unstable  for  sufficiently  narrow  distributions 
(see  record  521,12). 

The  program  also  calculates  "East  and  North"  statistics, 


covariance 


, M = - V E.  N.  - 


n i—  1 i 


IV  2 

standard  deviation  of  covariance,  °m  = ^ *EiNi^" 


T S 


standard  error  of  covariance 
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correlation  coefficient,  M 


M 


The  program  also  calculates  parameters  related  to  vector  quantities: 

the  scalar  amplitude  of  the  vector  mean,  = V E 2 + N2;  vector  variance 

V = — (G2  + O2) ; standard  deviation  = A/V  . 
v 2 E N ’ v 

Progressive  Vector  Diagram  (PROVEC) 

The  progressive  vector  displacements  are  plotted.  The  plot  begins 
with  an  asterisk  (*>  on  a day  boundary.  All  following  interval  boundaries 
are  indicated  by  a symbol,  which  changes  with  each  month  in  the  case  of 
longer  records . This  type  of  plot  accentuates  very  low  frequency  events 
at  the  expense  of  higher  frequency  oscillations. 

Variable  vs.  Time  Plot 

This  is  a diagram  of  any  variable  plotted  as  a function  of  TIME. 

The  plot  is  generated  from  the  1-hour  averaged  series.  This  type  of  plot 
is  complementary  to  the  progressive  vector  diagram  since  it  accentuates 
higher  frequency  events  such  as  inertial  and  tidal  oscillations. 

Spectra 

The  program  TIMSAN  (Time  Series  ANalysis)  uses  the  Fast  Fourier 
Transform  algorithm  of  Singleton  (1969)  and  is  restricted  to  data  segments 
of  length  N points,  where  N must  be  an  even  number  which  has  no  prime 
factor  larger  than  5,  and  must  be  less  than  8000  points;  data  series  longer 
than  this  must  be  broken  into  two  or  more  pieces. 

The  number  of  degrees  of  freedom  for  the  first  40  plotted  points  is 
given  by  v = a m s where  m is  the  number  of  adjacent  frequency  bands 
being  averaged  as  stated  in  the  label,  s is  the  number  of  independent 
data  pieces  being  averaged,  again  as  stated  in  the  label,  and  a should 
be  set  equal  to  2 for  temperature  spectra  and  for  Horizontal  Kinetic  Energy 
( HKE)  spectra  for  which  the  EAST  and  NORTH  components  seem  statistically 
independent.  In  the  absence  of  information  regarding  NORTH-EAST  cor- 
relation, one  should  use  a = 2 to  be  safe. 
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The  loq-loq  plot  is  further  averaged  during  plotting  so  that  more 


and  more  [>o i nts  art'  averaged  toqi'ther  as  frequency  increases.  This 
eliminates  the  bunchinq  toqetlier  of  jKtints  at  hiqh  f requenc i t's , increases 
the  degrees  of  freedom  of  the  high  frequency  estimates,  and  still  permits 
low-frequency  resolution.  The  averaging  algorithm  is  as  follows:  count- 
ing from  the  left  of  the  plot,  the  first  40  plotted  points  represent 
data  that  has  been  averaged  as  stated  in  the  label;  the  data  for  the  next 
15  plotted  points  has  been  averaged  over  twice  as  many  frequency  bands; 
the  next  6 over  five  times  as  many,  the  next  40  over  ten  times  as  many, 
the  next  15  over  twenty  times  as  many,  the  next  6 over  fifty  times  as 
many,  the  next  40  over  100  times  as  many  and  so  on.  In  this  way,  for 
example,  7900  data  points  with  no  averaging  indicated  in  the  label  would 
be  plotted  as  only  176  points,  and  the  last  14  estimates  would  be 
averaged  over  200  basic  frequency  bands.  The  m in  the  formula  V = a m s 
for  degrees  of  freedom  is,  in  this  example,  200  times  larger  at  the 
highest  frequencies  than  at  the  lowest  frequencies. 

For  V ^ 30,  the  confidence  limits  for  the  spectral  estimates  are 
given  approximately  by  (1  - 2/1V  ±z/2/9V) 5 where  Z = 1.28375  for 
80%  confidence  limits,  Z = 1.645  for  90%,  Z = 1.96  for  95%  and 
Z = 2.5757  for  99%.  In  the  example  above,  if  the  HKE  spectral  plot 
label  had  indicated  2 pieces  and  averaging  over  8 adjacent  frequency 
bands  then  v = 2 « 2 x 8 = 32  for  the  lowest  frequencies  (assuming 
NORTH  and  EAST  components  are  highly  correlated)  and  200  x 32  = 6400 
for  the  highest  frequencies.  The  95%  confidence  intervals  (i.e.,  95% 
of  the  time  one  would  expect  the  spectral  estimates  to  vary  no  more  than 
this  much)  would  be  (0.57,  1.55)  at  low  frequencies,  and  (0.97,  1.03) 
at  high  frequencies. 

For  v < 30,  one  must  obtain  confidence  intervals  from  Chi-Squared 
distribution  tables  in  standard  statistical  references. 

Temperature  spectra  from  850  current  meters  show  aliasing  at  high 
frequencies  because  of  the  sampling  scheme.  Speed  and  direction  are 
burst  sampled  within  each  recording  interval.  These  samples  are  vector 
averaged  to  attenuate  high  frequency  aliasing  in  the  current  measurement. 

There  is  only  one  5.19  second  temperature  sample  per  recording  interval, 
however.  The  temperature  spectrum  for  data  series  shows  clearly  the 
resulting  high  frequency  aliasing.  In  contrast,  the  VACM  averages  all 
variables  over  the  full  recording  interval,  virtually  eliminating  the  problem. 
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Gaussian  Filtered  Plots 


The  basic  time  series  is  filtered  using  a symmetrical  running 
Gaussian  filter  with  a half  width  of  24  hours.  The  filtering  is  sequen- 
tial and  the  resultant  time  series  is  48  hours  shorter  than  the  input 
time  series  (the  first  and  last  24  hours  are  lost) , and  the  sampling 
interval  is  the  same.  A simple  running  hat  filter  is  then  applied  to 
the  series  to  form  a series  with  one  data  point  per  day,  the  point  being 
the  average  for  each  day  measured  from  midnight  to  midnight. 

A)  Vector  Stick  Plots 

The  filtered  velocity  vectors  are  plotted  along  a time  scale. 
Unless  otherwise  indicated,  the  vector  orientation  is  such  that  north 
is  up.  The  length  of  the  "stick"  is  proportional  to  speed. 

B)  Line  Plots 

Pressure  and  temperature  are  filtered  as  described  above,  and  one 
point  per  day  is  plotted  along  a time  scale. 

C)  Sequential  Array  Plots 

Subsampled  filtered  velocity  vectors  or  temperature  values  are 
plotted  to  show  a spatial  array  at  a particular  depth. 
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